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ABSTRACT 

Gas motions in galaxy clusters play important roles in determining the properties of the intracluster medium 
(ICM) and constraining cosmological parameters using X-ray and Sunyaev-Zel’dovich effect observations of 
galaxy clusters. The upcoming ASTRO-H mission, equipped with high-resolution X-ray spectrometer, will 
make the first direct measurements of gas motions in galaxy clusters through measurements of Doppler shifting 
and broadening of emission lines. However, the physical interpretation of the data will be challenging due to 
the complex thermal and velocity structures of the ICM. In this work, we investigate how well we can measure 
bulk and turbulent gas motions in the ICM with ASTRO-H, by analyzing mock ASTRO-H simulations of 
galaxy clusters extracted from cosmological hydrodynamic simulations. We assess how photon counts, spectral 
fitting methods, multiphase ICM structure, deprojections, and region selection affect the measurements of gas 
motions. We show that while ASTRO-H is capable of recovering the underlying spherically averaged velocity 
profiles to within 20% with reasonable amount of photon counts (> 200) in the 6.7 keV Ee XXV line complex, 
there are considerable azimuthal variations in the ICM velocities, even in dynamically relaxed systems, which 
must be taken into account when interpreting data and developing observing strategies. Einally, we show that 
ASTRO-H should enable direct measurements of the hydrostatic mass bias with an accuracy of < 5%, by 
accounting for both rotational and random velocities from Doppler shifts and broadening of emission lines. 

Our results are broadly applicable for future X-ray missions, such as Athena+ and SMART-X. 

Subject headings: cosmology: theory - galaxies: clusters: general - methods: numerical - X-rays: galaxies: 
clusters 


1. INTRODUCTION 

In the hierarchical structure formation scenario, galaxy 
clusters are believed to grow through mergers and continuous 
mass accretion through the cosmic-web of large-scale struc¬ 
tures. These mergers and accretion events stir up the hot gas in 
the deep gravitational potential w ell of the clu ster, gen eratin g 
large-scale motions (Vazza et al. l2009Ll2011b Miniati l2015h . 
Non-thermal pressure provided by these internal gas motions 
exerts an extra non-thermal pressur e support to th e intra clus- 
ter mediu m (IC M) (Lau et al. 120091: Ba ttaglia et al. 1201 21 : Nel¬ 
son et al. l2014t Shi &. Komatsu l2014t Shi et al. 1201 5h .' intro¬ 
duces bias in the hydro static mass estimate s of galaxy c lusters 
(e.g., Rasia et al. I2006b Nagai et al. l2007bl : Nelson et al. l2012t 
Ettori et al. 1201 3h . and serves as one of the primary sources 
of systematic uncertaintie s in cl uster-based cos mologi cal con¬ 
straints (e.g., Allen et al. 120081 : Vikhlinin et al. 1200^ Planck 
Collaboration |20 1 5ll . However, the nature of gas motions in 
galaxy clusters are still poorly understood. While there are in¬ 
direct constraints on the level of gas motions thro ugh sp atial 
fluctuations in thermal pressure (Schuec ker et al. 120041) and 
X-ray surface brightness (Churazov et al. 120121) . direct obser¬ 
vational constraints from emission lines are stil l missing, and 
only u pper limits are av ailable (Ota et al. 120071 : Tamura et al. 
|2^ Sanders & Eabian lM^ . 

The upcoming ASTRO-H X-ray mission (Takahashi et al. 
I2014h . equipped with high energy resoluti on (^ 5 eV) Soft 
X-ray Spectrometer (SXS) (Mitsuda et al. 120141) . will make 
the first direct measurement of bulk and turbulent gas mo¬ 
tions, through shifting and broadening of the 6.7 keV Ee XXV 
Ka line. These measurements will provide new insights into 
the physics of energy feedback from active galactic nuclei 


(AGN), dynamics of cluster mer gers, a nd chemical enrich¬ 
ment of the ICM (Kitayama et al. 120141) . However, since the 
thermal and velocity structures of the ICM are rich and com¬ 
plex, interpretation of upcoming ASTRO-H measurements re¬ 
quires detailed cosmological simulations that include all rele¬ 
vant astrophysical processes. 

As the first measurements of ICM velocities are expected 
soon, there are several questions that must be addressed be¬ 
fore we can analyze and interpret the ASTRO-H data in a 
meaningful way. How well can ASTRO-H measure and sepa¬ 
rate bulk and turbulent motions in clusters? What are optimal 
observing and analysis strategies? Can we use ASTRO-H to 
measure the non-thermal pressure in the ICM and correct for 
the hydrostatic mass bias of galaxy clusters? 

The primary goal of this work is to address these questions 
by analyzing mock ASTRO-H observations of galaxy clusters 
extracted from high-resolution hydrodynamical cosmological 
simulations. These simulations capture internal motions in 
the ICM generated by mergers and accretion events in the hi¬ 
erarchical structure formation model. We create a realistic 
mock ASTRO-H data of simulated clusters, analyze them us¬ 
ing the data analysis pipeline used to analyze real X-ray data, 
and compare derived velocity measurements to “true” values 
measured directly in the simulations. We find that the spher¬ 
ically mass-weighted averaged ICM temperature and veloc¬ 
ity profiles can be measured with the accuracy of better than 
20%, which leads to good recovery of the hydrostatic mass 
bias to 5%, provided that there are at least 200 photons in the 
6.7 keV Ee XXV Ka line. However, our results also reveal 
considerable azimuthal variations in both bulk and turbulent 
velocities in the ICM (by up to factors of 2) at a given radius 
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even for dynamically relaxed looking clusters. We discuss 
implications of these azimuthal variations in correcting the 
hydrostatic mass bias using ASTRO-H cluster observations. 

This paper is organized as follows. Section|2]describes nu¬ 
merical simulations and mock ASTRO-H simulations used in 
this work. Results are presented in Section 3 We present 
discussion and conclusions in Section|4]and|5] respectively. 

2. SIMULATIONS 

2.1. Hydrodynamical Simulations 

The simulations presented in this work are performed using 
the Adaptive Refinement Tree (AR T) A-b ody-rgas-dy namic s 
code (Kravtsov ll999l : Kravtsov et al. l2002t Rudd et al. 120081) . 
which is an Eulerian code that uses adaptive refinement in 
space and ti me, and non-adaptive refinement in mass (Klypin 
et al. 1200 Ih to achieve the dynamic ranges necessary to re¬ 
solve the cores of halos formed in self-consistent cosmologi¬ 
cal simulations. The code is parallelized using Massage Pass¬ 
ing Interface (MPI) libraries and OpenMP directives. 

We analyzed a sample of high-resolution hydrodynamical 
simulat ions of galaxy clusters formation from Nagai et al. 
(l2007allR hereafter N07), performed with non-radiative gas 
physics. Likely targets for ASTRO-H would be nearby X- 
ray luminous massive clusters. Thus, we focus on two X- 
ray luminous clusters, CL104 and CLIOI, at z = 0, with the 
core-excised X-ray temperature of 7x = 7.7 and 8.7 keV, re¬ 
spectively. Each cluster is simulated using a 128^ uniform 
grid with 8 levels of refinement. Clusters are selected from 
120/i“' Mpc computational boxes, achieving peak spatial res¬ 
olution of sa 3 . 6 / 2 “' kpc, sufficient to resolve dense gas clumps 
in the ICM. The dark matter particle mass in the region sur¬ 
rounding the cluster is 9 x 10*/2“*Mq, while the rest of the 
simulation volume is followed with lower mass and spatial 
resolution. We refer readers to N07 for further details. CLlOl 
is a massive, dynamically active cluster, which has recently 
experienced violent mergers (at z ^ 0.1 and z ^ 0.25) and con¬ 
tains two major substructures near the core at z = 0. These two 
substructures have been identified by visual inspection and 
masked out before its analysis. CL104 is a similarly massive 
but dynamically relaxed cluster, with no significant mergers 
occurred in the last 8 Gyrs. 

In our previous work (Nagai et al. I2013h . we investigated 
effects of baryonic physics on the velocity structure of the 
ICM, by analyzing CL 104 and CLIOI performed with vary¬ 
ing gas physics. The first set is performed with non-radiative 
gas physics. The second set includes radiative cooling, star 
formation, metal enrichment, and stellar feedback (CSL). 
The third set includes CSL and energy feedback from AGN 
(CSLh-AGN). By comparing these three runs, we found that, 
while the baryonic physics can have significant impact on the 
properties of cluster cores, the velocity structure outside the 
cluster core (r > 0. ISrsoo) is relatively unaffected by the bary¬ 
onic physics. Since the present work concerns measurements 
of large-scale gas motions outside of cluster cores regions, we 
focus our analyses on the non-radiative simulations and dis¬ 
cuss their limitations and caveats in Section|4] 

2.2. Mock ASTRO-H Simulations 

We generate a mock ASTRO-H surface brightness map and 
spectrum for three orthogonal projection axes of each simu¬ 
lated cluster. The ASTRO-H pipeline consists of two main 
steps; (1) generation of the flux map from the simulation out¬ 
put; and (2) conversion of the flux map to a photon map by 



Fig. 1Mock ASTRO-H image of CL104 in the [0.3,10.0] keV band. 
9x9 boxes, each with a size of 3' X 3', were used for the analysis of projected 
velocity field. The grid (/, y) at the bottom-left comer is (1,1), and the grid 
label i j j indicates the cell located at the upward/rightward direction from the 
origin, respectively. 
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Fig. 2. — Mock ASTRO-H spectmm for the energy range [0.3,10.0] keV 
extracted from the grid (/,y) = (7,5) defined in Figure[T]of CL104 with total 
exposure time of 2 Ms. The inset shows the spectrum around the Fe XXV Ko 
line (redshifted to 6.27 keV) where the velocity constraints come from. The 
redline shows the best model BAPEC fit of the spectmm. 

convolving with the instrumental response of ASTRO-H. Lor 
the first step, we compute the emissivity ce as a function of 
energy E for each hydrod ynamic cell using the APEC plasma 
code (Srmthet al. 120011) with AtomDB yersion 2.0.2 (Los- 
ter et al. 120121) . The emissiyity ce = e£(p,r,Z,Zobs,vios) is a 
function of the gas density p, temperature T, the gas metal- 
licity Z = 0.3Zq, and line-of-sight yelocity vios of the cell. 
The obseryed redshift of the clusters is set to Zobs = 0.068, 
which corresponds to the redshift of A1795, a nearby relaxed 
cluster, which is one of the likely targets for mapping out the 
gas yelocity structure out to large radii with ASTRO-H. The 
energy range is £ S [0.1,10.0] keV with energy bin size of 
AL = 1 eV, and we include the effect of thermal broadening 
on the emission lines. 

Each flux map is computed by summing up the emis- 
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sivity of all cells contained within a sphere of radius 
of 5/i“* Mpc centered on the densest dark matter parti¬ 
cle. We then convolve our flux map with the ASTRO- 
H ARF (sxt-s_140505_ts02um_intall.arf) and 
RMF (ah_sxs_5ev_2 0130 80 6. rmf) response files from 
The energy resolution of the RMF file is 5 eV. We 
then draw photons for each location from the convolved flux 
map assuming a Poisson distribution. As we are interested 
in the intrinsic systematics, such as projection effects, rather 
than statistical uncertainties, we set the exposure time fexp = 
2 Ms per pointing to ensure enough photons in the spectra. We 
do not model background noise because it is subdominant to 
the strong Fe XXV Ka line where the gas velocity constraints 
come from. 


2.3. ASTRO-HData Analyses 

In order to assess how accurately ASTRO-H can measure 
temperature and velocity structure in the ICM, we analyze the 
mock spectra and study how the physical parameters depend 
on various analysis conditions, fitting method, energy band, 
region selection, and deprojection. Unbinned spectra (1 spec¬ 
tral bin = 1 eV) are analyzed using the Cash statistic, while the 
data are rebinned so that each spectral bin contains at least 25 
counts in the fitting. Considering the size of the SXS field 
of view, the spectral regions are defined by 3' x 3' grids (Fig¬ 
ure [D or annular rings with 3' width. The spectral fitting is 
performed by using XSPEC version 12.8. The spectral model 
consists of the BAPEC model and the Galactic absorp tion 
model (the wabs model; Morrison & McCammon ll983h and 
is convolved by the detector and telescope responses. The gas 
temperature, metal abundance, redshift, velocity dispersion, 
and the normalization are allowed to vary while the hydro¬ 
gen column density is fixed at Nh = 2 x 10 ^° cm“^, the value 
adopted when generating the mock simulatio ns 3 E or metal 
abundance, the tables in Anders & Grevesse (119891) is used. 
The velocity dispersion is defined by the Gaussian sigma for 
velocity broadening. The bulk velocity along the line of sight 
is calculated by c(z - Zciuster) where z is the fitted redshift and 
^cluster = 0.068. Eigure |2] shows a mock ASTRO-H spectrum 
extracted from the grid (/, j) = (7,5) in Eigure[T]of CL104 with 
total exposure time of 2 Msec. 

The deprojection of annular spectra to recover the three- 
dimensional source properties is performed by the pro jet 
model in XSPEC, assuming that the source emissivity is con¬ 
stant within the spherical shells whose radii correspond to the 
annuli used to extract the spectra. The quoted errors indicate 
Icr errors. 

To quantify the accuracies of velocity measurements from 
the mock data, we compare them against the “true" values 
measured directly from the simulations. Specifically, we com¬ 
pute the mass-weighted spherically averaged temperature and 
velocity profiles for each analyz ed clu ster following the pro¬ 
cedure described in Lau et al. (120091) . Briefly, we first de¬ 
fine the rest frame of the cluster to be the mass-weighted 
dark matter velocity within rsoo. We then compute the gas 
mass-weighted mean (v')mw and root-mean-square velocity 
in this rest-frame. We define the “turbulent” veloc¬ 
ity to be the velocity dispersion = y/We 

^ http://hea-WWW.harvard.edu/simx/ 

^ Tlie actua) va)ue of Nh toward A1795 is 1.2 X 10^® according to 
tlie LAB survey of Galactic HI. This difference in Nh has no consequence on 
our results and conclusions of the paper. 




o 



Fig. 3.— Temperature (top panel), bulk velocity (middle panel), and tur¬ 
bulent velocity (bottom panel) as a function of counts in the Fe XXV Ka line 
complex from fitting faked spectra of the isothermal gas with kT = 6.0 keV, 
Z = 0.3Zq in 5.0—10 keV with the BAPEC model. Different colors indicate 
different input turbulent velocity cturb = 0 (black), 100 (red), 200 (green), 300 
(blue), 400 (cyan), 500 (magenta), and 600 kms“* (orange). 


compute the “bulk” velocity as the projected mass-weighted 
velocity averaged along the line-of-sight. 

3. RESULTS 

3.1. Dependence on Photon Counts in the Fe XXVKa Line 

The primary constraints on the ICM velocity measurements 
with ASTRO-H come from the shifting and broadening of the 
bright 6.7 keV Ee XXV Ka line, therefore we first study the 
effect of photon counts in the Ee line on the parameter estima¬ 
tion. We simulated SXS spectra of isothermal gas assuming 
the BAPEC model with kT = 6.0 keV, Z = 0.3Zq, and input 
turbulent velocity of Uturb =0,100,... ,600 km s“* and zero 
bulk velocity. To study the effects of photon counts, we vary 
the normalization so that the Ee XXV Ka line complex con¬ 
tains 30 - 30000 counts. The simulated spectra in 5-10 keV 
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Fig. 4.— Left panels: results of the BAPEC model fitting to the 5.0-10.0 keV spectra using the Cash statistic for CL104. Middle panels: ratios of BAPEC 
model parameters from spectral fitting using the Cash and statistic in the 5.0-10 keV band for CL104. Right panels: ratio of BAPEC model parameters from 
spectral fitting between 5.0-10 keV and 0.3 — 10 keV. The temperature (top panels), the bulk velocity (middle panels), and the turbulent velocity (bottom panels) 
are plotted along the i-axis of the map shown in Figure[T] 


were fitted by the BAPEC model. Figure [3] shows the best-fit 
parameters and their Icr statistical errors using Cash statistic. 
It shows that > 200 counts in the Fe XXV Ka line complex are 
required to measure the turbulent velocity accurately through 
line broadening for the assumed range of velocities. The er¬ 
ror range in Figure [3 was derived by first simulating a large 
number of spectra and fitted the input model to each and com¬ 
plied a histogram of the best-fit values from each trial and es¬ 
timate the mean and the root-mean-square value. To achieve 
the accuracy better than 20%, > 200 counts are needed in the 
6.7 keV Fe XXV Ka line complex when the turbulent velocity 
Cturb < 200 kms“'. We recommend > 200 counts in 6.7 keV 
Fe XXV Ka line complex for robust measurements of (Ttuib- 
Note that the requirement for calibration accuracy of SXS 
energy-scale is 2 eV with a goal of 1 eV (Mitsuda et al. l2014h . 
which corresponds to the line shift of 90 km s“' (or 45 km s“') 
at 6.7 keV. Thus, in very bright regions like cluster cores, the 
measurement is limited by the instrumental uncertainty rather 
than the statistical error. 

3.2. Fitting Statistics: Cash vi. 

Two statistics. Cash (Cash Il979h and are commonly 
used in X-ray spectral fitting. We compare how these two 


statistics differ in the parameter estimation for the relaxed 
cluster CF104. Because velocity measurement requires suffi¬ 
cient photon counts (> 200 counts) in the Fe XXV Ka line, the 
spectral regions are limited to the central 5x5 grids (approx¬ 
imately r < r 25 oo). The left panels of Figure |4] shows results 
of the BAPEC model fitting to the 5.0-10.0 keV spectra us¬ 
ing the Cash statistic and the middle panels shows comparison 
of resultant parameters between Cash and statistics. Each 
spectrum contains at least 200 counts in the Fe XXV Ka line. 
The differences in the fitted temperature, bulk and turbulent 
velocities between the two different methods are typically less 
than 20%, which is reasonable as the Cash statistic approaches 
the x^ statistic in the limit of large photon counts. Since two 
methods produce consistent results, the unbinned spectra are 
analyzed by utilizing the Cash statistic in the following anal¬ 
ysis. Cash statistic is preferred for low photon counts, the 
spectrum does not need to be rebinned as in the case of x^ 
statistic, thus preserving the spectral energy resolution crucial 
for measuring gas velocities from line broadening. 

3.3. Spectral Energy bands and Multiphase ICM 

The temperature of the ICM may affect the accuracies 
of turbulent velocity measurements from line broadening. 
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Fig. 5. — Comparison of bulk (top) and turbulent (bottom) velocities de¬ 
rived in the fiducial case with the true gas temperature and the fixed tempera¬ 
ture case with Tx = 7.7 keV, for CL104. Different color points correspond to 
different j location in the map shown in Figure[T] 

The line broadening originates from both thermal and tur¬ 
bulent motions and given by a combination of Aftth = 
EoikT and Afsturb = £o<7turbC~*, where Eq is the rest- 

frame energy of the line emission, and m is the ion mass. 
Since the temperature kT is determined primarily by the shape 
of continuum spectrum of the hot ICM, it is ideal to use a 
wider energy range for accurate parameter estimation. If a 
certain amount of cool gas exists along the line of sight, the 
soft emission may affect the temperature determination even 
though the emissivity of Fe XXV Ka line at low temperature is 
low. We thus compare analysis in two different energy ranges, 
5.0-10 keV and 0.3-10 keV, where the latter includes more 
emission from cool gas. 

As seen from the right panels in Figure |4] for the relaxed 
cluster CL 104, the fitted parameters do not depend on the 
choice of the energy range for most regions except for the 
grid {i,j) = (7,7), where the temperature derived from 5.0- 
10 keV is ^ 25% higher than that derived from 0.3-10 keV. 
This is due to the presence of an infalling cooler substructure 
which is visible in the north-west of the cluster image, shown 
in Figure [T] 

The recovered turbulent velocity does not show clear de¬ 
pendence on the energy range. This is expected as most 
of the constraints on the turbulent velocity come from the 
broadening of the Fe XXV Ka line present in both energy 
ranges. Since the velocity constraints from line broadening 
is inversely proportional to the mass of the ion species, ve¬ 
locity measurements based on the Doppler broadening of the 
Fe XXV Ka ion is least affected by thermal broadening com¬ 
pared to emission lines of species of lighter elements (e.g.. 



Fig. 6.— Top panels: spectral regions for the deprojection analysis for 
CL104: 6 annuli (left) and 4 boxes in 8 azimuthal directions (right). Bot¬ 
tom panels: the same regions are shown for CLIOI. The maps show the 
jc-projections for t he two clusters. The location of merger shock is indicated 
by the arrow (see ^3.4.3F 

Inogamov & Sunvaev [2003h . On the other hand, the soft band 
(0.3-5.0 keV) contains emission lines from lighter elements 
which are more susceptible to thermal broadening, and leads 
to biases in turbulent velocity estimates. 

We further examine the effects of multiphase ICM by fix¬ 
ing the gas temperature to 7.7 keV, while keeping the gas den¬ 
sity, abundance, exposure time and fitting method unchanged. 
The velocities in this case should be unbiased from the effects 
of the multi-temperature ICM. We compare the velocity es¬ 
timates between the fixed temperature case and the fiducial, 
multi-temperature case in Figure |5] While the bulk veloc¬ 
ity measurements of the two cases are consistent with each 
other within statistical errorbars, the turbulent velocity is un¬ 
derestimated by ^ 30% in the fiducial case compared to the 
fixed temperature case. This is expected because the turbu¬ 
lent velocity (determined from broadening of the Fe XXV Ka 
line) is more affected by higher temperature gas through ther¬ 
mal broadening, while the bulk velocity (determined from line 
shift) is not. 

3.4. Velocity Profiles 
3.4.1. Deprojection analysis 

One of the key science goals of ASTRO-H is to measure 
non-thermal pressure profile due to bulk and turbulent gas 
motions and their impact on the hydrostatic mass. The non- 
thermal pressure profile is determined from spherically aver¬ 
aged density and velocity profiles. To study how well we can 
recover three-dimensional temperature and velocity profiles 
from the ASTRO-H spectra, we extract 0.3- 10 keV spectra 
from annular rings with 3' width (top left panel in Figure |6l) 
and fitted them with the BAPEC model. 

Figure |2] show the projected and deprojected quantities for 
the X, y, z projections for a relaxed cluster CL 104. Since 
CL 104 is a relatively relaxed, spherically symmetric cluster, 
the temperature and turbulent and bulk velocities from the 
three orthogonal projections agree with each other to within 
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Fig. 7.— Top panels: Deprojected temperature (left panels), bulk velocity (middle panels), and turbulent velocity (right panels) of the relaxed cluster CL104 
viewed along three orthogonal projections: x (black), y (red), z (green) projections (open circles). The projected quantities are shown with the crosses for 
comparison. Bottom panels: Deprojected temperature, bulk velocity, and turbulent velocity for the 8 azimuthal directions (crosses) for the x projection. The Cash 
statistic was used to fit the 0.3 — 10 keV spectra. The solid lines indicate mass-weighted, spherically averaged values computed directly from the simulation, and 
the dotted line indicates r 25 oo of the cluster. 
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Fig. 8.— Effects of sample size. Results of deprojection analysis for spectra extracted from 1 arm (W), 2 arms (W-hN, W-hE, W-i-S), 4 arms (W-i-N-i-E-i-S), 
and 8 ai'ms are plotted. Deprojected temperature (top panels), bulk velocity (middle panels), and turbulent velocity (bottom panels) of the relaxed cluster CL 104 
viewed along the x projection (crosses). The Cash statistic was used to fit the 0.3-10 keV spectra. The circles indicate deprojected values. The solid lines indicate 
mass-weighted, spherically averaged values computed directly from simulation, and the dotted line indicates r 25 oo of the cluster. 


about 10%. Furthermore, we find that the deprojected values 
are in agreement with the mass-weighted, spherically aver¬ 
aged profiles measured directly from the simulation to within 
< 20%. Note, however, that the recovered turbulent velocity 
tends to be underestimated by 10% in some annuli, particu¬ 
larly in the outer regions (r > ^2500) where the photon counts 
in the Fe XXV Ka line drops below 200. 

3.4.2. Azimuthal Variations 

Since the Field-Of-View (FOV) of SXS is small, we need to 
perform mosaic observations with multiple pointings to cover 
the spatially extended regions at large cluster-centric radii. 
However, since accurate velocity measurements require fairly 
deep exposures (> 100 ks at r « r 25 oo). it would be costly 
to perform mosaic observations. In practice, observations at 
the large cluster-centric radii will be limited a few directions. 


Since the cluster gas is not perfectly spherically symmetric 
and its velocity distribution is generally anisotropic, we ex¬ 
pect that the velocity profile measured along any given direc¬ 
tion is not representative of the azimuthally averaged velocity 
profile in an annular ring. 

In order to investigate the level of the azimuthal variations 
in the recovered velocity profiles, we extract spectra from four 
grids at r < ^2500 in eight azimuthal directions: 0°, 45°, 90°, 
135°, 180°, 225°, 270°, and 315°from the horizontal, corre¬ 
sponding to W, NW, N, NE, E, SE, S, and SW (top right panel 
in Eigure|6]l. The deprojected temperature and velocity pro¬ 
files of the 8 azimuthal directions in the x projection is plotted 
in Eigure|7]for CL104 in the left and right panels respectively. 

Even in the relaxed cluster CL104, we find significant scat¬ 
ter in the fitted temperature, bulk and turbulent velocities 
among different azimuthal directions. The azimuthal varia- 
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Fig. 9.— Same as Figure|2]for the merging cluster CLIOI. 


tion in the projected temperature increases with radius from 
10% at r = 0.2r25oo to 30% at r = ^ 2500 - The temperature varia¬ 
tion in the outer bin is driven by the infalling cold substructure 
in the NW direction. The same substructure is also partly re¬ 
sponsible for the large differences in the bulk and turbulent 
velocities. 

These results suggest that the velocity profile measured 
along any single arm is not representative of the spheri¬ 
cally averaged profile. Specifically, the deprojected turbulent 
velocities along any single direction systematically under¬ 
estimate the spherically averaged turbulent velocity measured 
directly from the simulation (shown in the bottom-right panel 
of Figure |2ll, because the simple spherical analysis fails to 
capture the azimuthal variations in the bulk velocities (shown 
the bottom-middle panel of Figure |2l). Such variations in the 
bulk velocities are expected, because mergers and filamentary 
accretions occur anisotropically along the axes of mergers and 
filaments, respectively. 

In general, decomposition of gas motions into turbulent and 
bulk components, inferred from broadening and shifting of X- 
ray emission lines respectively, depends on the size and geom¬ 
etry of the regions where the X-ray spectra were extracted (see 
the Appendix for more discussion). To study the dependence 
of velocity recovery on the size of the regions, we co-added 
spectra at the same radii in different pointing directions (e.g., 
W and N) and performed the deprojection analysis. Figure [8] 
shows the resulting recovered temperature, bulk, and turbu¬ 
lent velocity when the spectra in 2,4 or 8 arms are combined. 
Compared to the results of a single arm, the turbulent/bulk ve¬ 
locity converges to the spherically averaged values when we 
combine spectra from more arms. This suggests that sam¬ 
pling more multiple directions leads more robust estimates of 
the velocity profiles. Having 8 arms leads to better recovery 
of the velocity profiles to about 20%. 


Next, we repeat the same deprojection and azimuthal anal¬ 
yses on the merging cluster CLlOl. The bottom left and right 
panels of Figure |6] show the regions selected for these analy¬ 
ses, respectively. Figure |9] shows its deprojected profiles for 
the temperature, bulk and turbulent velocities in annular bins 
and in the 8 azimuthal directions, compared with the simula¬ 
tion values. In general, the temperature and turbulent veloc¬ 
ities from the annular analysis recover the true values quite 
well, although the scatter among different orthogonal projec¬ 
tions are greater than that of CL104. The recovery of the bulk 
velocities is generally accurate to within 20% at r < r 25 ooc, 
but it is underestimated at larger radii where there is a merger 
shock (the approximate location is indicated in Fig. |6]l. The 
azimuthal variations in both temperature and velocities are 
also larger for the unrelaxed cluster CLIOI. 

3.5. Hydrostatic Mass Bias 

One of the primary sources of astrophysical uncertainties 
in cluster cosmology lies in hydrostatic mass bias due to non- 
thermal pressure support provided by internal gas motions in 
the ICM. Assuming spherical symmetry, the total cluster mass 
within some radius r can be expressed as the following; 

M(< r) = MHSE+47i-and+4^rot+4^accel- (1) 

where Mhse is the hydrostatic mass, M^nd is the mass term 
associated with random (or turbulent) gas motions, Mmi is the 
term associated with the rotational gas motions, and Maccei 
is the term associa ted wi th acceleration (or deceleration) of 
gas (see Lau et al. 1201 3L for details). Note that other mass 
terms are negligible (< 1 %) co mpare d to M^nd, Mmt (Lau et al. 
120091) andMaccei (Nelson et al. 120141) . The hydrostatic mass is 
given by 


-rkT f d\ogp cHogT^ 
Cpnip \91ogr d\ogr ) ’ 


3.4.3. Dependence on Dynamical States 


(2) 
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Fig. 10.— Recovery of the hydrostatic mass bias in the relaxed cluster 
CL 104 (top panel) and the merging cluster CLIOI (bottom panel). Upper 
panels: The black line represents the true mass profile measured directly from 
simulation, the red points represent the hydrostatic mass Mhse measured 
from simulation, and the magenta points represent the recovered mass profile, 
Mrec =^HSE+^rand+^rot, using the velocity measurements from ASTRO-H 
mocks. Lower panels: Profiles of the hydrostatic mass bias, MnsE/^true - 1 
(red points) and the recovered mass, jV/rec/^true -1 (magenta points). Differ¬ 
ent symbols indicate the recovered mass Mrec obtained from the ASTRO-H 
mock maps viewed along three orthogonal projections of each cluster. The 
vertical dotted lines indicate r 25 ooc for oach cluster. 


where T and p are the mass-weighted gas temperature and gas 
density, respectively. Note that the hydrostatic mass depends 
on the gradients of gas temperature and gas density, which 
may not be well determined with ASTRO-H owing to its rel¬ 
atively poor spatial resolution of 1'3. The mass terms due to 
velocity Miand and Mrot, however, are less sensitive to resolu¬ 
tion because both the average bulk and turbulent velocities do 
not change significantly with cluster-centric radius, and thus 


Fig. 11 .— Histograms of the mass bias of the recovered mass, Mrec/Mtrue — 
1, at r = 0.5 Mpc, computed from the deprojected bulk and turbulent velocity 
measurements along individual azimuthal arms for the relaxed cluster CL 104 
(top panel) and the merging cluster CLlOl (bottom panel). The magenta 
lines represent the recovered mass computed from combined turbulent and 
bulk motion measurements, while the the blue and green lines represent the 
mass computed using only random and bulk motions, respectively. 


can be determined rather well with ASTRO-H; 
_-r^dpa^ r 
Gp dr '■ 

r 2 

Mml = - p , (Vrot) , 

o 


(3) 

(4) 


where tr^ are the mass-weighted velocity dispersion 

for the radial and tangential component, and (viot) is the mass- 
weighted mean rotational velocity. The computation of these 
mass terms requires three-dimensional velocity fields that are 
not directly accessible to ASTRO-H observations. Only ob- 
servationally accessible quantities are line-of-sight measure¬ 
ments of turbulent velocity (ctlos) and bulk velocity (vlos) 
measured using broadening and shifting of X-ray emission 
lines, respectively. To estimate M^nd from observational data. 
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we assume the gas velocity dispersion is isotropic, mean¬ 
ing that the measured line-of-sight turbulent velocity is equal 
to the isotropic one-dimensional velocity dispersion: CTlos = 
(T[d = = ct?/ 2- Note that hydrodynamical simulations sug¬ 

gest that, while gas velocity dispersion is isotropic around our 
radius of interest (^2500 ^ r < r^Qo) on average, it can be fairly 
anisotropic for individual clusters with 0 = 1 -afjilaj) « 
±0.5 (e.g., see Figure 1 in Lau et al. I2009h . leading to scatter 
in Mrand- For Mrot, we estimate the rotational velocity (vrot) 
as the difference in the line-of-sight bulk velocity (vlos) be¬ 
tween two opposite regions; e.g., (vrot) ~ |vlos,nw-vlos,se| 
where vlos.nw and vlos.se are line-of-sight bulk velocities in 
the NW and SE directions. Note that this estimate of the ro¬ 
tational velocity assumes that the axis of rotation lies in the 
plane of the sky, which is not necessarily true. 

The upper panels of Figure[T0]show the mass profiles of the 
two clusters: CL104 (top panel) and CLIOI (bottom panel). 
The black lines represent the true mass profile measured di¬ 
rectly from simulation, the red lines represent the hydrostatic 
mass Mhse measured from simulation, and the magenta points 
represent the recovered mass Mrec = MnsE+A/rand+A/i-ot, with 
the latter two terms computed with velocity profile measure¬ 
ments from ASTRO-H mocks. Different symbols indicate 
the recovered mass Mrec obtained from the ASTRO-H mock 
maps viewed along three orthogonal projections of each clus¬ 
ter. The lower panels of Figure [TO] show the profiles of the 
HSE mass bias and the recovered mass, relative to the true 
cluster mass, for the same clusters. The HSE mass bias is 
shown in red, and the bias of the recovered masses Mrec = 
fHHSE+fHrand+Afrot is shown in magenta. 

In the lower panels, we demonstrate that ASTRO-H clus¬ 
ter observations can be used to recover the total cluster mass 
at r > r 25 oo at the level of about 5% for dynamically re¬ 
laxed clusters, provided that we measure and correct for both 
fHrand (random motion) and Mrot (rotational motion) terms (see 
the appendix for more discussion). For example, for the re¬ 
laxed cluster CF104, the remaining bias in the corrected mass, 
Mrec/Mtrue-1, Varies from-1-1% to-4% in the outer three bins, 
compared to -3% to 7% for the HSE mass. For the unre¬ 
laxed cluster CFlOl, the improvement in mass measurements 
is more evident; i.e., the remaining bias in the corrected mass, 
Mrec/Mtrue - F improves from -9% to -1-3% in the inner bins. 
However, the mass recovery is not as good in the outer bins, 
because the gas is not in s teady- state. 

As shown in Section 13.4.21 the velocity measurements 
along individual azimuthal arms are not representative of 
spherically averaged velocity profiles, which could introduce 
biases in the recovered mass. Here, we demonstrate that this 
mass bias due to azimuthal variations in velocity can be miti¬ 
gated provided that we account for both bulk and random and 
turbulent gas motions from line shifts and broadening, respec¬ 
tively. In Figure [TTl we show the histograms of the recovered 
mass computed from the deprojected bulk and turbulent ve¬ 
locities in each azimuthal arm, compared to true cluster mass 
at r = 0.5 Mpc. Note that we have made the assumption that 
the rotational velocity is equal to the deprojected bulk veloc¬ 
ity. For the relaxed cluster CF104, the recovered mass for a 
single arm computed from only turbulent (bulk) motions un¬ 
derestimates the true mass by -4 (-1)%. The mass recovery 
improves when we include both random and bulk motions, 
leading to smaller biases: -3% to -1-2%. For the merging clus¬ 
ter CFlOl, the mass recovery from a single arm also improves 
if we consider both bulk and turbulent motions. 


4. DISCUSSION 

There are several effects that are missing from the present 
mock ASTRO-H analysis. We briefly discuss their effects on 
our results. 

Point Spread Function: A systematic effect may arise from 
scattering of photons by a relatively wide Point Spread Func¬ 
tion (PSF) of the SXT onboard ASTR O-H. The half power 
diameter is 1.3' (Takahashi et al. 120141) . Strongly peaked X- 
ray emission at the center of relaxed cluster may contaminate 
the measurement of turbulent velocity at fainter, outer regions. 
Abundance peak in cool core clusters will result in stronger 
contamination. In the case of A1795, the fraction of scat¬ 
tered photons in a SXS field of view placed at the scale ra¬ 
dius r 25 oo from adjacent regions is estimated to be 25% based 
on SXS simulations assuming the azimuthally-averaged PSF 
model. This effect can be corrected by simultaneously analyz¬ 
ing a series of pointings along the radius from the center and 
solving intrinsic distributions of the line shifts and width (see, 
ASTRO-H White Paper o n Clus ters of Galaxies and Related 
Science by Kitayama et al. l2014h . 

Baryonic physics: The results presented in this paper are 
based on non-radiative simulations that do not include bary¬ 
onic effects such as radiative cooling, star formation, and 
AGN feedback. These baryonic physics can alter the velocity 
structure of the ICM especially in c luster cores, where their 
effects are significant (Nagai et al. I2013h . Radiative cool¬ 
ing and star formation produces the multiphase ICM stric¬ 
ture through thermal instability (e.g., Voit & Donahue l2015h . 
but the cooling flow is regulated by peri odic bu rsts of out¬ 
flows from the central AGN (e.g., Fabian 120121 : Zhuravleva 
et al. I2014h . generating the complex ICM structure in this re¬ 
gion. However, the velocity structure outside the cluster core 
(r > 0.15r^nn ) is rel atively unaffected by the baryonic physics 
(Nagai et al. l2013h . Furthermore, the ICM temperature has 
shown to be more homogeneous outside the clus ter cor e with 
the inclusion of baryon dissipation (Rasia et al. I2014h . sug¬ 
gesting that the systematic uncertainties due to temperature 
inhomogeneities should be smaller. 

Resonant Scattering: We did not include the effects of res- 
onant scattering, e. g., in the Fe XXV line (Zhuravleva et al. 
I 2 OI lb Shang & Oh l2013h . Resonant scattering increases the 
optical depth and suppresses the strength of the emission line. 
The effect of resonant scattering can be important in most cen¬ 
tral region of the cluster where the density, and therefore the 
optical depth to the 6.7 keV Fe XXV Ka line is high. Out¬ 
side the central region, however the effect is expected to be 
small. In the presence of turbulent motions, resonant scat¬ 
tering is suppressed by Doppler shifting the line out of reso¬ 
nance. Adding the effect of resonant scattering should there¬ 
fore improve the velocity recovery of ASTRO-H as it provides 
additional constr aints f rom the shape of the resonant line (e.g., 
Zhuravleva et al. l2013h . 

5. CONCLUSIONS 

Gas motions play an important role in determining the prop¬ 
erties of the hot X-ray emitting ICM and constraining cosmo¬ 
logical parameters using X-ray and SZ effect observations of 
galaxy clusters. Despite their importance, we know very little 
about the nature of gas motions in the ICM. The upcoming 
ASTRO-H mission will provide the first direct measurements 
of gas motions in galaxy clusters through Doppler shifting and 
broadening of X-ray emission lines. In this work, we investi- 
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gate how well ASTRO-H can measure the bulk and turbulent 
gas motions in the ICM by analyzing mock ASTRO-H maps 
and spectra of simulated clusters extracted from hydrodynam- 
ical cosmological simulations. We summarize our main find¬ 
ings below: 

1. ASTRO-H measurements of the bulk and random ve¬ 
locities are accurate to < 20%, when the number of 
photons exceeds 200 in the 6.7 keV Fe XXV Ka line 
complex. This enables ASTRO-H to make velocity ve¬ 
locity measurements out to r < r 25 oo along any given 
arm of a nearby rich galaxy clusters with a reasonable 
exposure time (< 100 ks). However, ASTRO-H veloc¬ 
ity measurements become increasingly more expensive, 
requiring deep exposure (> 1 Ms) in the low surface 
brightness region in the outskirts of galaxy clusters. 

2. Velocity measurements using only the hard (5.0- 
10 keV) band of the X-ray spectrum achieve similar ac¬ 
curacies as those based on the full (0.3- 10 keV) band 
as both bands contain the 6.7 keV Fe XXV Ka line. 
However, using the soft (0.3-5.0 keV) band spectrum 
alone introduces a significant bias, as the emission lines 
in this band come from low mass ions that are more 
susceptible to thermal broadening, making the effect of 
Doppler broadening more difficult to discern. 

3. Spectra extracted from annular bins give a good recov¬ 
ery of the thermodynamic and velocity profiles of the 
ICM. Specifically, deprojected temperature and veloc¬ 
ity profiles from spectra extracted from the annular bins 
agree well with the “true” mass-weighted, spherically 
averaged profiles computed directly from the simula¬ 
tions to better than 20%. 

4. We find significant azimuthal variations in the recov¬ 
ered turbulent velocities, by up to a factor of 2, even 
in dynamically relaxed clusters. Our finding suggests 
that the velocity structure of the ICM is richer and more 
complex than its thermodynamic counterpart. Measure¬ 
ment of the turbulent velocity profile in any given az¬ 
imuthal direction (or arm) is, therefore, not representa¬ 
tive of the underlying mass-weighted, spherically aver¬ 
aged profile, which is required to correct for the hydro- 
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Static mass bias. Since the field of view of ASTRO-H 
SXS is small, mosaic observations are required to sam¬ 
ple the annular regions at the large cluster-centric radii 
of nearby clusters, which could be expensive and infea¬ 
sible. Our results suggest that combining spectra from 
eight arms leads to better recovery of the velocity pro¬ 
files to about 20%. 

5. We demonstrate that ASTRO-H observations can be 
used to recover the total cluster mass at r < r 25 oo at 
the level of 5% for dynamically relaxed clusters, pro¬ 
vided that we measure and correct for mass biases as¬ 
sociated with both random and rotational motions. The 
remaining uncertainties arise from the assumptions on 
the isotropy of random gas motions and projection ef¬ 
fects on recovered bulk velocities. Note that accurate 
measurements of the hydrostatic mass require high- 
angular resolution X-ray instruments, such as Chandra 
and XMM-Newton, and ASTRO-H (with limited angu¬ 
lar resolution) provides the hydrostatic mass bias asso¬ 
ciated with bulk and random gas motions. 

Although our analyses have focused on ASTRO-H, some of 
these results are broadly applicable to future X-ray missions, 
such as Athena-^ and SMART-]^ Future missions will im¬ 
prove sensitivities and angular resolutions, produce spatially- 
resolved velocity structure in the ICM, and enable measure- 
ments of the ICM velocity power spectrum (Zhuravleva et al. 
I2012h . However, future measurements will be increasingly 
limited by systematic uncertainties (associated e.g., with mul¬ 
tiphase ICM and dynamical state of clusters) rather than sta¬ 
tistical uncertainties. Future work must address how to sepa¬ 
rate bulk and turbulent motions, how well do different codes 
agree on the gas velocity field, and what are the main drivers 
of bulk and turbulence motions. 
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APPENDIX 

DECOMPOSITION OF BULK AND RANDOM GAS MOTIONS AND HYDROSTATIC MASS BIAS 

In Lau et al. (1201 3h . we showed that there are two methods for measuring hydrostatic mass bias from gas motions, and that they 
are mathematically equivalent. Here we give a brief review of the two methods. The first method, called the “Summation method”, 
is based on Gauss’s Law and Euler’s Equation, where the cluster mass within some volume V with surface dV. Assuming steady- 
state of gas (for simplicity), the mass estimate obtained using the “Summation method” is given by 

(— + {yV)y)-dS, (Al) 

47rG Jdv\ P J 

where v is the gas velocity on the infinitesimal surface element dS, and p and P are the corresponding gas density and pressure. 
We can recover the true mas s of the cluster if we measure v, p and P of individual gas element and sum them up over the entire 
surface using Equation (lAlb . In reality, it is not possible to do so since observations are limited by finite spatial resolution of 
the instrument and we cannot measure the physical quantities of each individual gas element. In p ractice, we measure spatially 
averaged values above some resolution scale. In this spirit, we can spatially average Equation (lAll) over the surface dV to arrive 
at a cluster mass estimate that is based on spatially averages of gas velocity, density and pressure. The mass obtained using this 
“Averaging method” is given by 

where the angular bracket (•••) denotes spatial averaging, and denotes spatial averaging weighted by the gas density p 

over a spherical shell with radius r. This mass estimate is the same as , since mass is an extensive quantity and th erefo re it 
is not affected by any spatial averaging. Averaging of the nonlinear term (v • V)v introduces to a new term in Equation (IA2b that 
contains the density-weighted velocity dispersion tensor a}j = ((v', - (vi)p)ivj - {vj)p))p. This term represents random motions of 
the gas elements internal to the resolution scale, which manifests as Doppler broadening in the emission lines. The bulk motion 
of the gas at th e res olution scale is represented by (v)p, which can be measured as Doppler shifts in the emission lines. Thus, 
using Equation (IA2b . we can recover cluster mass by accounting for the internal random motions and bulk motions from Doppler 
line broadening and shift, respectively. In the limit of large averaging scale, the motions within the volume are uncorrelated, and 
the random motion term dominates. At very small scale where we can resolve motions of individual gas elements, the motions 
consist solely of the bulk motion of the gas element and the random motion term is zero. While the breakdown of the motions into 
“bulk” and “random” components depends on the averaging scale, the recovered mass from these gas motions is independent 
of the averaging scale. This can be understood by considering the fact that the Summation method is a limiting case of the 
Averaging method where the averaging scale approaches zero. Since both methods give the same mass estimates, the recovered 
mass from the Averaging method for any arbitrary averaging scale should give the same cluster mass as the Summation method. 
This means that as long as we measure both shifting and broadening of the emission line to account for both bulk and random 
gas motions, we can recover the cluster mass independent of the resolution scales of instruments. 







